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Abstract 



Irrotational relativistic vortex configurations in uniform subsonic motion 
with respect to a surrounding perfect fluid are analysed for the purpose of 
, application to superfluid layers in neutron stars. Asymptotic solutions are 

I found by asymptotically expanding the flow equation at large distances from 

' the vortex core and then by solving it order by order. The asymptotic effec- 

■ five tension and energy density that are needed for an averaged macroscopic 

description are thus obtained as functions of the vortex velocity, the vortex 
circulation, the asymptotic chemical potential and of parameters depending 
on the equation of state. 



o 

^ ■ I. INTRODUCTION 

> ' 

^ , The aim of this work is to obtain large scale average values of physically important quan- 

^ I titles such as energy and tension for stationary longitudinally invariant "vortex" solutions 
for a relativistic irrotational perfect fluid in an asymptotically uniform background in a flat 
spacetime. The present analysis is restricted to the case in which only a single vortex is 
present. This analysis constitutes an essential first step towards the large scale treatment of 
cases in which, as discussed in the concluding section, there is an extended array of vortices. 

The main application we have in mind is the interior layers of a rotating neutron star, in 
which the bulk of the matter can be represented to a good approximation as a relativistic 
superfluid at zero temperature, whose local behaviour on a microscopic scale is that of an 
irrotational perfect fluid. This means that the vorticity resulting from the stellar rotation 
will be concentrated in discrete vortex cores which may be subject to external drag forces 
that cause them to move relative to the ambient superfluid. In order to evaluate the effect 
of such relative movement, the present study generalises results of an earlier investigation 
that was restricted to the axisymmetric case in which the vortex is at rest with respect 
to the uniform background, so that an exact analytic description was possible even for the 
inner regions of the vortex. In that previous work, the fluid under investigation was a 
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superfluid at non zero temperature, which can be described by the relativistic extension 
of Landau's two fluid model [^]. Thanks to the numerous symmetries of the configuration, 
namely invariance under time translation, invariance under longitudinal translation (along 
the axis of the vortex) and axisymmetry, it was possible to obtain exact analytical results 
for a "cool" superfluid (Bose condensate plus a perfect gas of phonons 0). Here, we still 
assume invariance under time translation and longitudinal translation but axisymmetry is 
lost because of the motion of the vortex with respect to the surrounding fluid. 

This loss of axisymmetry makes it difficult to obtain exact analytic results except in the 
extreme limit of a "stiff" fluid - in which the speed of sound is equal to that of light - which is 
relevant for various cosmological contexts, notably the case in which the superfluid is formed 
as an axionic fleld condensate whose vortices constitute what are known as "global" cosmic 
strings 0. The present work is intended for the rather different astrophysical context of 
neutron star interiors - in which the relevant sound speeds are expected to be considerably 
(albeit not incomparably) smaller than the speed of light. In these circumstances it is not 
easy to provide an exact analytic treatment of the inner regions of a moving vortex. However 
the present work shows that it it is possible to provide a very satisfactory description of the 
outer regions by an asymptotic expansion. This description is quite sufficient for the purpose 
of calculating the large scale averages (of energy, tension et cetera) that are needed in practice 
for the usual astrophysical applications. 

It will be convenient to work with cylindrical coordinates with respect to which the 
spacetime metric gf^u (p, ci = 0, 1, 2, 3) is given by 

Qp^dx^dx" = -c^dt^ + df + dr^ + r^dO^ , (1.1) 

where c is the speed of light, so that the postulated stationarity and longitudinal invariance 
will be expressible as the condition that the relevant physical quantities should be indepen- 
dent of t and i. It is to be remarked that the neglect of general relativistic curvature effects 
in this work is justifled in so far as most neutron star (and even cosmological) applications 
are concerned by the consideration that the gravitational curvature radii involved will be 
extremely large compared with the lengthscales - of which the most important are the inter- 
vortex separation distances - that are relevant for the local effects considered here. It will be 
supposed that the latter distances are themselves very large compared with the lengthscales 
characterising the central vortex core, whose details may consequently be neglected without 
signiflcant loss of accuracy. (Neglect of the core would not of course be justiflable if we were 
concerned with cosmic strings of "local" as opposed to "global" type, nor of the analogous 
phenomen of the quantised magnetic flux tubes in neutron star matter.) 

The present analysis will not allow for thermal effects, whose treatment would require 
the use of the generalised Landau type two constituent superfluid model 0, IQ, but whose 
consequences are not very important in typical neutron star applications for which the 
zero temperature limit approximation is sufficient. Our attention will be restricted here to 
superfluid models of the simplest "barotropic" kind for which a complete set of independent 
physical variables is constituted just by the components of the 4-momentum covector fi^ 
whose dynamics is governed by the irrotationality condition 

/ia = V,</?, (1.2) 

where is a gauge dependent scalar - interpretable as a phase angle in the underlying 
bosonic condensate - with period 2tt in units such that the Dirac Plank constant h is set 
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to unity. The complete system of the equations of motion depends on the specification of 
the appropriate equation of state, giving the fluid pressure P as a function of the square 
of the effective mass or "relativistic chemical potential" fi - defined as the magnitude of 
the momentum covector - which determines the corresponding dilatonic amplitude scalar $ 
according to the prescription 

2 dP 

= ^1-2 ' = -/i>a = -(/'.aVV. (1.3) 

This in turn determines the corresponding particle current vector n°", which is given by 

nP = $V , (1.4) 
while the associated stress energy momentum density tensor will be given by 

T"" = $VV + Pg""" ■ (1-5) 

In conjunction with (1.2), all that is needed to complete the system of dynamical equations 
of the superfiuid (thereby automatically ensuring the conservation of the stress momentum 
energy tensor) is the equation expressing the conservation just of the current (1.4). This 
will evidently be presentable as a non-linear wave equation in the form 

V.($V» = . (1.6) 

The special "stiff' (Zel'dovich type) case is characterised by the special condition that <l>^ 
simply be constant (which arises for an equation of state of the form P oc /i^ — m^, where m 
is a fixed mass per particle), so that ( p..6|) reduces to a linear equation of the familiar Dalem- 
bertian form, for which the bicharacteristic speed Cg of propagation of small perturbations 
is evidently equal to the speed of light c. For the more general class of equations of state 
covered by the present treatment, the corresponding bicharacteristic sound speed Cg will be 
given by the formula 

For an equation of state of "relativistic polytropic" type P oc /i^"^^ where w is a fixed 
"polytropic index" number, one obtains the fixed value (Cg/c)^ = 1/ n. 

In order to determine the asymptotically contributions to the modifications of energy and 
tension due to the presence of a vortex, the only information about the uniform background 
state that was found to be needed in the axisymmetric non-moving case consisted of 
the asymptotic limit values, /x^ and say, of the squared chemical potential /i^ and the 
squared amplitude <I>^, together with the corresponding limit value of the squared sound 
speed Cg^ which it will be convenient to express in terms of a dimensionless non-negative 
parameter 5 defined as a measure of deviation from "stiffness" by the formula 



^4((e/cJ-l)^(^|,)L. (1.) 



In the case of a simple relativistic polytrope it will be given by 25 = w — 1, which means 
that it will be equal to unity, 5 = 1, in the case w = 3 for which the equation of state is that 
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of the standard model for a gas of ultrarclativistic particles, which has been commonly used 
as a crude first approximation for the treatment of neutron star matter. In the case of a 
relatively moving mortex, it will be shown here that this information is not quite sufficient, 
but that it is also necessary to know one more parameter, characterising the next order 
of differentiation of the equation of state. This extra parameter, 5'' say, is conveniently 
specifiable in the form 

so as to be interpretable as a dimensionless measure of deviation from the relativistic poly- 
tropic case. Although it is conceivable that it might be negative, one would expect in practice 
that for reahstic equations of state this parameter 5'' would usually be positive . 

In addition to this minimal information about the uniform background state, the only 
other parameter needed to characterise the vortex in the non- moving axisymmetric case is 
the corresponding circulation integral k which (since we have set h = 1) must be a multiple 
of 27r, or equivalently the corresponding angular momentum A say per unit mass, taking 
the latter to be the asymptotic value ii^ of the relativistic chemical potential /i, provided it 
is taken for granted that the value of the net current input / per unit length is zero. These 
conserved quantities will be expressible as averages over a circle of radius r (whose value 
can be chosen arbitrarily without affecting the result) by the formulae 

A = , K=(fd(p = 2^ (1.10) 

and 

/ = 2^rrr?I^, (1.11) 

where indicates angular averaging over the chosen circle of radius r, and a prime indi- 
cates differentation with respect to the angle 9 round the circle, while i^o-, with components 
(0,0,1,0) in the system (1.1), is the outgoing unit normal to the circle. Finally, in the case 
of a moving vortex, the specification of A , and if necessary of / , must evidently be supple- 
mented by the specification of the velocity, f3 say, of the relative motion in order for the 
characterisation of the vortex to be complete. It is conceivable in principle that / might be 
given a non-zero value by an artificially contrived injection process in a laboratory experi- 
ment, but in the natural context of neutron star vortices it is reasonable to assume that - 
as will be postulated in the last sections of the present work - the central core input should 
vanish: / = 0. 

The main results of the present work are the calculation of the asymptotic form of the 
tension, i.e. the integrated longitudinal stress, T say, and the corresponding energy per unit 
length U say, over a circular cross section of radius r through the vortex, as functions of the 
background parameters fi^, 6 and 5^ and of the parameters A and P that specify the 
amplitude and velocity of the vortex itself, subject to the assumption that the injection rate 
/ is zero. As in the familiar non- moving case these quantities are found to have a logarithmic 
radial dependence. 



(1.9) 
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II. THE ASYMPTOTIC EXPANSION AND ITS LOWEST ORDER 



Our objective is to evaluate large scale averages of energy and tension for stationary, 
longitudinally invariant, asymptotically uniform solutions for an irrotational barotropic rel- 
ativistic perfect fluid in flat space. It is not possible to solve explicitly the superfluid dy- 
namical equation ( |1.6| ) except in the special "stiff' case where this equation becomes linear. 
Since we are ultimately interested by the asymptotic effect of the vortex and not the de- 
tails of the flow near the core of the vortex, we shall try to solve this equation by using an 
asymptotic expansion at large r. Prior to this, it is convenient to reexpress ( |1.6|) in the form 

2V.V> + ((C/C3) 2 - l)/i-2(v,^2)y.^ ^ Q _ (2.1) 

At this stage, one can invoke the time translation symmetry and the longitudinal sym- 
metry, to which one can associate two Killing vectors k"^ and P generating respectively time 
and longitudinal space translations and corresponding to the ignorable coordinates t and i. 
These two Killing vectors provide us with two constants of motion 

ky^ = -E, ly^ = c, (2.2) 

where E is interpretable as the effective energy per particle, and C is interpretable as an 
effective longitudinal momentum per particle. Consequently, the potential (f has a linear 
dependence on time t and longitudinal distance i in cylindrical coordinates. Since there is 
no loss of generality in eliminating the £ dependence altogether, i.e. setting £ = 0, by a 
longitudinal boost, the potential will be conveniently expressible in the form 



^(4-0 (2.3) 



where u is a function only of r and 6. 

In order to analyse that asymptotic form of the flow at large r, we now postulate that u 
has an expansion of the form 



u = vr -\- w 



Inr z /1\ , , 

Inr + x+y + - + 0- , 2.4 



where 1/, x, y, and z are functions of 6 only. In order to proceed, we need the first and 

second partial derivatives of u with respect to r and with respect to 9. To the relevant 
asymptotic order, the former will be expressible as 

Ur = v + --y— + '-^ + o[ — ), 2.5 

w In r 2z — 3y / 1 \ , , 

= --2+'^y— + + ' (2-6) 
while the latter will be expressible as 

/ /, / /Inr z' /1\ 

ug = ^'r + w'lnr + x' + y' + - + - , 2.7 



v"r + w" In r + x" + y"— + — + o(-), (2.8) 



5 



using a prime to denote straight differentiation with respect to 9. 
According to (|1.3|) , the effective mass /i will be given by 

cV' = ^(c^ - u,' - ^,u,') . (2.9) 

This can be evaluated with the help of the asymptotic forms of the first partial derivatives 
( [^.5D and ( |2.7| ) and we obtain 

cV^ 2 / 2 , /2n o ' '^"^^ (vw + vV) ;2(ln!£ 

— r- = c — (l/ +1/ j — 2l/ w 2 w — 

II //sin'" liv z — v' z' — V — x'^ — w*^ /1\ 
+2(.y - v'y' - ^'.')— + ^ -fl + o(-) . (2.10) 

As well as the expression for /i^ itself, we shall use its partial derivatives with respect to r 
and 9 for the evaluation to the required order of both the quadratic first order differential 
contribution 

(/i^),,V> = ^{{^^\ Ur + \{^l')e ue) , (2.11) 



and the second order Laplacian contribution 

E 11 
V^V"^ = -^(urr + -Ur + ^Uee) , (2.12) 

in the dynamical equation ( |2.1| ). 

Independently of the particular form of the equation of state, it can be seen that the 
vanishing of the coefficient of the leading term on the left of ( |2.1j ) , namely the term of order 
1/r, will give the requirement 

v" + v=0, (2.13) 



whose first integral gives 



v' + v'' = V\ (2.14) 



where is a constant of integration that will be interpretable as the flow velocity of the 
uniform background. There will be no loss of generality in arranging for the direction of 
this flow to be aligned with the axis 9 = 0, which means that the local solution for v will 
be given explicitly by 

v = Vcos9. (2.15) 

Still independently of the particular form of the equation of state function given by (p.. 71), 
it can also be seen that the vanishing of the coefficient of the next leading term on the left 
of ( p.l|) , namely the term of order In r/r^, will give the further requirement 

w' = (2.16) 

as the only possibility compatible with global regularity, which simply means that (like (3) 
w must be a constant. We shall restrict our attention in subsequent sections to the cases in 
which w is restricted to be zero, which is the condition for there to be no source or sink in 
the vortex, i.e. for which there is no net inflow or outflow, but for the purposes of the next 
section we shall provisionally continue to consider configurations of the most general kind 
for which the constant w is unrestricted. 
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III. CONDITIONS FOR EQUILIBRIUM AT INTERMEDIATE ORDER. 



In order to proceed beyond the lowest order, we must take account of the particular form 
of the function given by ( |1.7| ), which will have an asymptotic expansion that can be seen 
from (|2.10|) to take the form 

l((c/cJ^-l)=. + 24(^^^ + o(i). (3.1) 

using the usual abbreviation 7^ = 1/(1 —/S^) with f3 = V /c, where S and 5" are dimensionless 
constants defined in (|1.8|) and ( [L.9|) in terms of asymptotic values of respectively second and 
third derivatives of the equation of state function. 

On substitution of the leading order dynamical conditions (|2.14|) and ( p.l6| ), the expan- 
sion formula ( p. 101 ) simplifies to 



(i) ■ (3-2) 



cV^ ^wv + x'v' ,\nr 1/1/2 + x'2 + 2yi/ + 2(zi/')' 

It is now apparent from ( p.l| ) that the vanishing of the coefficient of the leading survivimg 
term the left of ( p.l|) , namely the term of order l/r^, will give the requirement 

- 251/ '2)x" + 4<5v/x' + 2./.<5(i/2 - i/'2) = . (3.3) 

This provides a decoupled second order differential equation for x, for which one immediately 
obtains a first integral expressible as 

( = A—^fY^ , (3.4) 

where A is an arbitrary constant of integration which is interpretable as proportional to the 
corresponding value of the total circulation k round the vortex, in terms of which it is given 
by the relation 



K = f dip , (3.5) 



27iE ' / 
in which a is a predetermined constant given by 

which, assuming the asymptotic fiow velocity V = cj3 to be subsonic, will necessarily be less 
than unity: 

V^<c^ ^ a^<l. (3.7) 

Subject to this subsonicity condition, which will be taken for granted in all that follows, 
the first order differential equation (p.4|) can in its turn be integrated straightforewardly to 
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give the explicit solution for x which is found to be simply expressible in terms of a modified 
angle variable by 



X = A'j + wln\/l — a'^ sin^ 6 , ip = arctanjVl — tan 6'} (3-8) 

The absence of a second arbitrary constant of integration from this expression does not 
imply any loss of generality, because the absence of radial dependence in the corresponding 
term in ( |2.4| ) means that it is only the derivative of x but not its absolute value that is 
physically relevant. For the same reason (unlike the other angular dependence functions 
V, w, y, z which must all be strictly continouous) the function x can be discontinuous- 
and indeed it necessarily will be so somewhere unless the circulation n vanishes. It can 
be observed that while the physically meaningful derivative x' is a continuous function of 
6, the expression (|3.8|) gives a value of the potential x itself that has jump discontinuities 
at ^ = ±71/2. (Insertion and adjustment of an additative constant of integration on the 
right of (|3.8| ) could be used to make gauge changes of ip by which the discontinities could 
be displaced elsewhere, or by which one - though not both - of them could be removed 
altogether.) 

Again using (|3.1|) , we can go on to evaluate the coefficient of the next leading term in 
( p.l| ), namely the one of order In r/r^, whose vanishing can be seen to be expressible as the 
condition 

- 26v'^)y" + 86vv'y' + + 26{v'^ - 2v'))y = . (3.9) 

We thus again obtain a decoupled second order equation, this time for y, for which it can 
be seen that there is a first integral given by 

2 2 

(l - 25^/2) v'y' + (l + 25^v'^) vy = VB , (3.10) 

where 6 is a constant of integration. The complete solution is hence finally obtainable in 
the explicit form 

B cos 6 + C sin 6 

y = — ' (3-11) 

1 — a"^ sm u 

where C is another constant of integration. The regularity of this solution is guaranteed 
by the subsonicity condition (|3.71 ) which ensures that the denominator does not vanish 
anywhere. It is to be noted that unlike the solution ( |3.8| ) for x, the solution ( |3.11 ) for y is 
not affected by the presence of a non vanishing source coefficient w. 

Although the values of the integration constants C and 6 appear at this stage to be 
arbitrary, it transpires that they are in fact predetermined by the condition of existence of a 
globally regular solution for the still unknown function z . The latter is governed by a higher 
order equilibrium equation which we must now work out, even though it turns out that 
the integrals we shall need for evaluating the asymptotically dominant contribution to the 
energy and tension can be evaluated without any specific knowledge of the actual functional 
form of the higher order coefficient z itself. 
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IV. CONDITIONS FOR EQUILIBRIUM AT THE RELEVANT HIGHER ORDER. 

In view of the automatic cancellation expressed by ( [4.1[ ), the dominant surviving con- 
tribution to the integrals with which we are concerned will come from the next term in 
( p.2|) , which is only of order l/r^. In order to obtain the information required for evaluat- 
ing this contribution, we now consider the next leading term in ( |2.1| ) which (since we have 
already dealt with all the terms up to In r/r^) is the one of order The vanishing of the 
corresponding coefficient, as evaluated using (p.l|), is found to be expressible as the condition 

- 25.y' + 85./z' + + 26iv- - 2.2))z - 2(^ + 5(.- - 3.^))y 

1 / /2\ /2 



-46vv'y' - 2(1 + 25 + 25Vc-2/2) v'x'x" + 4(5 + 25^^/2) v> 

2 2 
= 2vx"w + 45^^(v/x" + (/2 _ 3^2)x')/u. - 4(5 + 5^\{v^ - v'))vw^ . (4.1) 



in which the right hand side will vanish when the source coefficient w is set to zero. 

For the sake of simplicity, let us from now on restrict our attention to the physically 
important source free case, meaning that we shall impose the postulate 

1^ = . (4.2) 



In these circumstances the right hand side of ([4.1|) will drop out, leaving an equation for z 
that can be rewritten, with the unknown terms grouped on the left, in the form 

2 2 / 2 

^-25v'%z' -{^~25{y^ + v^))v'^ = ((l + 25 + 25^\v'^)vv'x'^ + A5v\'y 



^2 J \^2 ' J \S (? 

2 2 

,2 _ ^/2 _ 2y 2|^25 + 5^^/2))x'2 + 2(^ - 5(V 2 + 2/2))^y . (4.3) 

The y dependent term and the term proportional to x'^ at the end of the preceeding ex- 
pression can be evaluated in the form of a derivative. We can thus obtain a first integrated 
version of the equation ( [4.3|) for z in the form 

(1 - sin^ ^) cos z' + (1 - a2(l + cos^ ^)) sin z = , (4.4) 

in terms of a known source function JF, from which the fully integrated solution will be 
obtainable the form 



/"cos 6* 



(4.5) 



\ — o? sin^ Q ' 

where the function f is calculable by direct quadrature as a solution of the differential relation 

^' = ^ • (4-6) 
cos^ y 



It can be seen that the inhomogeneous source term on the right of ( |4.4| ) will be given by an 
expression of the form 
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T = v + Q + n , 



(4.7) 



in which the first term is an automatically periodic symmetric contribution expressible by 

P = K + CflnVl-a'sin^e / , . — ^) 4.8 

^ 1 - q;2 sm^ Q ' 



where /< is a new constant of integration, while the second term is an antisymmetric contri- 
bution, with a manifestly periodic form expressible by 

l^/\2sin2g / l + /3V(2,52-^ + ,5^) _ (1 - a2)(l + 2(5 + 2^V^2 sin'g) \ 
1? 1^ 2(1 -a2 sin^^) 2(1 -«2 ^^^^ ) 

l-a^Q + cos^^) ,^ , 

+Bsin2g ^ ^ . 4.9 

2 1 — a"' sm 



What matters most for our present purpose is the third term in ([4.7|) , which is an antisym- 
metric remainder with a not automatically periodic form given by 

/^V(25-/3V(252-5-5^)) 



n = /(^^i\.yi2 ^^ + ev/r^^. (4.10) 

Unlike the other terms, the second contribution 7^ will be subject, for generic parameter 
values, to a mismatch between its value at departure from the axis of motion, when 

'?/' = 0, and its value at return to this axis, when = 27r, since it is evident from ( 4.10 ) that 
we shall have 

[^]o =2^i ^2(1 _ ^2)1/2 + eVl-«2J. (4.11) 

The only physically admissible solutions are those for which z and hence also T satisfies 
the condition of being a regular periodic function of the angle Q. Since this condition is 
automatically satisfied by the other contributions V and Q, it must therefore also be satisfied 
by the third term 'R. in ( [4.7|) , which means that the mismatch ( 4.11 ) must be made to vanish. 
It can thus be seen that the integration constant B cannot be chosen independently but that 
- as a necessary condition for global regularity - it must be specified in terms of the constant 
A (which according to (|3.5|) is proportional to the circulation k) by the relation 



B^Wl)-P^-W'-^-^>)). (4.12) 

^ c ^ (1 — a^) 

This not only eliminates the mismatch but actually gets rid of the remainder term altogether, 
so that we have 

7^ = 0, T = V^Q, (4.13) 

and it can also be seen to have the important secondary effect of removing the singularity 
that would otherwise have arisen from the the contribution to z via ( ^.5| ) and ( [4.6| ) from the 
antisymmetric term Q which, by substitution of ( |4.12|) , is rather miraculously simplified so 
as to be given by 
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(4.14) 



cos2^ c2 V (l-a2sin2^)2 2(1 -«2)(1 -a2gin2^^) y- 

It is evident that the corresponding contribution to the quadrature for f as given by ( [4 .61) 
will thus be satisfactorily regular and periodic. It is to be noted however that in order for 
z as given by ( [4.5|) to be appropriately regular and periodic, while it is of course necessary 
that f should also be periodic, it is nevertheless not necessary that f should be regular in 
the strictest sense, since it it is permissible for it to have a simple pole where cos 9 vanishes. 
Such a pole will in general arise from the contribution to ( [4.5|) from P, which can be seen 
to be expressible in the form 

^ ^an0(/< + Cln7l-a2sin2 0)-C^^=y . (4.15) 



cos2^ V ' ^ ' 



a 



In addition to the admissible pole term proportional to tan 9, the differentiand on the right of 
( [4.15| ) also contains an inadmissibbly non-periodic term proportional to the modified angle 
coordinate ip whenever the coefficient C is non zero. It can thus be seen that as well as fixing 
the constant of integration B by the condition ( [4.12| ) the requirement of global regularity 



also fixes the other constant of integration introduced in ( 3.11| ) by the simpler condition 

C = . (4.16) 
The solution for z that is thus finally obtained from ( [4.5| ) and ( |4.6| ) will be given explicitly 

by 



VA^cos9 / /3V(2(5' + ^ + ^^) (2^ - /?'7'(2^' -5- 5^)) In Vl - sin^ 9 ' 
^2 1^ a^il-a"^ sin2^)2 ^ {I - a^){l - sm^ 9) ^ 

K sin 9 + L cos 9 



1 — sin^ 9 
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(4.17) 



where L is yet another constant of integration, which, like K, presumably depends on the 
nature of the inner core. However what matters for the purpose of the following sections 
is not the particular functional form of the higher order angular coefficient z, but just the 
conditions ( |4.12| ) and ( [4.16| ) that are necessary and sufficient for its global regularity, and 



that specify the specific form of the intermediate order angular coefficient y which is what 
we really need to know. 



V. ASYMPTOTIC AVERAGES. 

Our purpose is to obtain asymptotic averages of quantities Q that will of the kind dis- 
cussed in the appendix of the preceeding work in the sense that - as in the prototype 
case of the pressure P - they can be postulated to have a smooth functional dependence 
on the chemical potential, and will thus have a deviation from the corresponding uniform 
asymptotic limit value Q^o that will have a form given by 



dQ_ 
dij? 



1 <PQ 

2 (c//i2)2 



2 \2 



(5.1) 
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The relevant deviation of the square of the chemical potential (or effective mass) variable n 
will itself be given, according to ( p.2|) , by 



This differs substantially from the corresponding expression in our previous analysis in 
which both v and y vanished (thereby also incidentally eliminating the term involving z) so 
that the leading terms of order 1/r and of order In r/r^ were absent, leaving just the term of 
order with coefficient simply given by x'^. Its effect is to introduce new leading terms 
of order 1/r and of order In r /r^ in the resulting expression for the asymptotic deviation of 
the quantity Q, which will be given by 

Q-Q =-t^2dQ (yx' ^ o(^.^.ylnr ^ x^^ + 2.y + 2(/zy N ^ ^^^4 d'Q W-'f , / 1 \ 

(5.3) 

The long range leading order terms in (|5.3|) will however cancel out, leaving a non zero 
contribution only from the short range term of order 1/r^ when we take the relevant average 
Q, as evaluated over a circular section of radius r through the vortex. Such an average will 
be given by an expression of the form 

Q = ^ f^''dr = '-^Q^ + ^ £^rdr (5.4) 

where is the average over an inner core of radius and where, for any given radius r, 
the quantity Q is the corresponding angular average, as given by 

"Q=^iQde. (5.5) 

Since the coefficient of term of order 1/r is antisymmetric with respect to 6 its angular 
average will evidently vanish. The angular averages of the term of order In r/r^ and of 
the term involving z will also vanish because they have the form of derivatives of periodic 
fuctions. The angular average of the deviation in which we are interested will thus be given 
just by 

n~"7r (,2dQ 7'v '2^^ o, 4 d^Q f^a\ 



Since the preceeding outcome is only of order 1/r^, it can be seen that the required 
average will have an asymptotic behaviour of exactly the same kind as was obtained in the 
previous analysis p|] of the non-moving case, i.e. it will be described by an expression of the 
form 



^ ^ Q \nr /lnr\ 

Q-Q^ = ^ — + o{—) , (5.7) 
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with a constant coefficient Q which will have the same dimensionality as the quantity Q itself, 
where /4 is a constant normalisation factor, interpretable as the average angular momentum 
per unit mass, that is specified independently of Q by the formula 

(5.8) 



which means that it is identifiable with the integration constant introduced in (|3.5| ). The 
explicit formula for the asymptotic deviation coefficient Q that is defined in this way will 
however be significantly more complicated in general than in the non-moving case, for which 
it was found to be given simply hy Q = —2{fi^dQ/dfi'^)\^ 0. The generalisation that is 
needed for the case of a moving vortex, i.e. one with a non-zero value of (3, can be read out 
from (62) as 

with dimensionless coefficients 3^ and Z given by 

To evaluate these coefficients, we use the explicit expressions obtained from ( |3.4] ) and from 
( CT ) and (CT) - We get 

^= 2(l-a2)i/2 ' (5-11) 

and 

(i-vT^)(a^-a^ + 7V + 2/3V^^) 

^ 2(l-a2)i/2 + a\\-a^) ' ^ ' ^ 

The latter expression is rather unweildy even in the "polytropic" case characterised by = 0, 
but in the stiff limit characterised by 5 = 5^ = it simply gives 3^ = 1, independently of the 
velocity /?. 

The applications with which we are principally concerned are the pressure P for which 
( pl?D gives 

P = -c'fil^liy-2ZS) , (5.13) 

and the squared amplitude $^ itself for which ( |5.3|) gives 

$^2 = _2$2 ^ys - 2Z{6^ -6- 6^)) . (5.14) 

These quantities are all that is needed to evaluate the average deviations of the longitudinal 
components 

Ui = (5.15) 
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of the conserved current and 



Tij = ^"^HiHj + Pg 



(5.16) 



of the stress energy momentum tensor, where the Latin indices run over the values 0,1 in 
the coordinate system ( |1 . 1| ) , with respect to which the gradient components 



(5.17) 



take the constant values 



fi, = -E , /i, = , (5.18) 

so that we simply obtain 

Hi = ^^fii (5.19) 

and 

fij = S'^fi.fij + Pgij. (5.20) 

The general expression ( |5.7| ) is not limited to scalar functions of yU^ but can be also 
applied to other quantities. This will however require more work because the intermedi- 
ate formulas such as ( p.6| ) are no longer applicable. Interesting examples of this kind are 
the transverse current components and the transverse-transverse components of the stress- 
momentum tensor. Although the resulting asymptotic coefficients are mathematically well 
defined, the physical interpretation of these quantities is less obvious than the components 
treated above so that we have preferred to relegate these results in an appendix. 



VI. CONCLUSIONS 



The present work has given asymptotic solutions for a relativistic vortex moving with 
respect to a surrounding irrotational perfect fluid (or, equivalently, a zero temperature su- 
perfluid). We have obtained the effective energy density, which can be expressed in the 
form 



eff\ 



4vr2 - 



Inr /Inr 



(6.1) 



and the effective tension, which can be given by 



U, 



eff\ 



_$2 

47,2 oo 



Inr 



dnr 



+ (6.2) 



where 3^ and Z are given by ( [5.11| ) and ( |5.12| ). 

The present work constitutes an intermediate step that is needed for the construction 
of a realistic modelization of the relativistic dynamics of a neutron star interior. Like a 
superfluid in a rotating bucket, the neutron superfluid in a rotating neutron star minimizes 
its free energy by forming an array of quantized vortices that, on a scale larger than the 
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intervortex separation length, simulates a rigid body rotation. Because of the huge number 
of vortices in a neutron star, it is not practicable to keep track of all individual vortices, so 
a macroscopic description requires the use of a model involving a vorticity 2-form 

= 2V[p/i^], (6.3) 

that is not supposed to be the small scale local vorticity field, (which vanishes outside the 
microscopic vortex cores) but that is to be interpreted as the large scale average over a 
neighbourhood extending across a many vortices. 

The previous analysis of a relativistic vortex configuration in a background without 
relative motion provided the basis for the specification of a particularly simple model p| 
within the general category that is needed for such a description. Such models are deriv- 
able from a Lagrangian type master function A{n'', Wpo-} that depends on the three scalar 
quantities that can be built from the (macroscopically averaged) current density n'' and the 
(macroscopically averaged) vorticity Wpo-, namely the magnitude n of the particle current 
vector n'' itself, the scalar magnitude w of the vorticity covector Wpa and the magnitude C 
of the associated Joukowsky lift force density vector C,p as defined by 

n2 = -n^n, , = ]^Wp„w'"' , C' = CpC^ (6-4) 

where the Joukowsky vector is defined as 

Cp = Wp„7f . (6.5) 

This vector is interpretable as representing the volume density of force that would be exerted 
on the vortices as an expression of the Magnus effect, by the relative flux of the fluid according 
to the simple formula originally derived by Joukowsky for flow past a long aerofoil. The 
earlier analysis was sufficient to unambiguously determine the appropriate form for the 
equation of state for the function A{n, w, (^} only in the limit for which C, vanishes. 

As convenient ansatz for use as a provisional approximation, it has been suggested |^ 
that this special limit form should be extrapolated by assuming that there is no dependence 
on C, at all, even when C, is non-zero, a supposition that is mathematically justifiable for a 
fiuid obeying the "stiff" Zel'dovich type equation of state, but not in general. To obtain a 
more accurate treatment for an arbitrarily compressible fluid it will ultimately be necessary 
to use a more sophisticated ansatz taking account of the effect of the relative flow that will 
be present when C, is non-zero. The present analysis of the effect of relative flow on an 
individual vortex is an indispensible flrststep towards the achievement of what is required. 
What still remains to be done is to extrapolate the present analysis to the case in which 
there is not just one vortex but an extended array of parallely aligned vortices. In the 
absence of relative flow at large distance, i.e. in the case corresponding to vanishing such 
an extrapolation was almost trivial, at least in the large separation limit that is relevant, 
because of the axial symmetry of the individual vortex solution. However for the non 
axisymmetric conflgurations considered in the present work the required extrapolation will 
not be so straightforward. 

It is to be remarked that the effect analysed in the present work is not the only kind 
that needs to be examined in greater detail in order to improve the accuracy of the recently 
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proposed model ^ for the treatment of relative flow past the vortices in neutron stars. 
Whereas the modification of the stress - momentum - energy tensor considered here is likely 
to be very small, a potentially more important correction arises from the friction between 
the vortices and the "normal" matter that will be present, not just due to thermal effects 
(which in typical neutron stars will be very small) but also, even in the zero temperature 
limit, due to the fact that the matter will not be entirely constituted just by neutrons, but 
will also include a small but significant fraction in the form of protons, together with a gas of 
charge neutralising electrons whose scattering from the vortices can provide a dynamically 
important friction drag []T|. 
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VII. APPENDIX: ASYMPTOTIC AVERAGES OF TRANSVERSE CURRENT 

COMPONENTS. 

In this appendix, we compute the asymptotic formula ( |5.7| ) for quantities which are not 
scalar functions of /i^ and therefore for which the method used in the main body does not 
apply. The first important example is that of the tranverse components with respect to 
Minkowskian coordinates characterised by 

x° = t , x'^ = i , = rcos6 , x^ = rsin^ (7.1) 

of the conserved current Since these components have no dependence on the longitudinal 
coordinates {i = 0, 1) the current conservation law automatically reduces from four to 
two dimensional form, 

V^n" = Van" = , (7.2) 

so that we obtain 

Vfe(xV) = , (7.3) 

using early Latin letters for the transverse indices, a,b = 2, 3. For any 2-dimensional trans- 
verse section S bounded by a curve s with outgoing unit normal having components i^a, 
Green's theorem thus provides us with the relation 

/ n"(iE = Ix^'n^ui^ds . (7.4) 

is Js 

In the case of a circular section with boundary characterised by 

x"- = ru"- , = cos9 , = sm9 (7.5) 

so that ds = rdO, the preceeding formula can be used to see that the corresponding average 
- defined as in ( |5.4|) - will be expressible just in terms of an angular average - as defined in 



(5^) - over the section boundary in the form 
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^=2z/Vfen'' . (7.6) 
In the same way, for the stress energy momentum tensor we obtain 

V^T% = VaT% = 0, (7.7) 

and hence 



Using 
we obtain 



T; = 2z/VbT^^ . (7.8) 



E 

Ua = , f^aT^"" = " "r- (7.9) 



.*;(-24^:^-(,+24(./)V)!^^)+o(i). (7,10) 



Since it is evident from symmetry considerations that the angular average x w z/" must 
vanish, we are left with a relation of the standard form ( p.7|) , namely 



n"- — n"' = n' 



with 



^/A^lnr /lnr\ 



2 2 

^ = -^o.^ 2 (y + 25^(y/)V)z/» , (7.12) 

in which the components of the relevant angular average are calculable from ( p.l5|) and 

( p. lip . Using the formula ( ^.12| ) for B, and noting that C vanishes by ( [4.16| ), one finally 



obtains 



2(25- /?V(252- 5-5^)) / 



a 



n^ = - -[1 h", ny = 0. (7.13) 

Having arrived at this result, it behoves us to remark however that the physical signif- 
icance of the transverse current component average thus obtained is rather less obvious 
than that of the corresponding longitudinal average obtained in the previous section, 
which is directly interpretable as being proportional to the current flux through the circular 
section across the cylindrical vortex cell under consideration. The corresponding transverse 
parts of the current flux associated with cylindrical vortex cells in a honeycomb lattice are 
also of obvious physical interest, but they are not directly definable as averages of the kind 
considered above, so their evaluation will require a rather different kind of calculation that 
will be left for a future article. 

The same question of physical significance arises for the averages of the transverse com- 
ponents of the stress energy momentum tensor, for which an analysis of the same kind as 
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that given in the preceeding paragraph also leads to an equation of the standard form, 
namely 

n-T:,-f-/-l'^.o{^). (7.14) 

In the particular case of the the cross components T"* this is obvious because they are 
expressible in the form 

T'** = /iV = nV (7.15) 

in which the longitudinal momentum components /i* are constant. It is thus evident that the 
corresponding averages - which have a direct physical interpretation in terms of longitudinal 
fluxes of transverse momentum - will simply be proportional to those of the current, being 
given by 

~ ^ _2(2S - m - S - / _ 1 - VT^y^, p ig) 

For the derivation of the corresponding formulae for the purely trans verse components T\ 
- which do not have such a straightforward physical interpretation - a little more work is 



required since (7^) no longer applies. One must then consider the product of the asymptotic 
expansion of given by ([5.41), with that of the given by ( |2.5D and (PT7|), in order to 
obtain an asymptotic expansion at the required order for the part $2^^^^ of the components 
Tab, the pressure part being already known via (|5.13| ). One obtains finally 

+ + (1 _ VT^^){2 + /3V -a^- 2(/3V^V«^)) 

(7.17) 

1^ = 0, (7.18) 



00 00 



yl — 



^ 2 2^2 (l-v^l^)(«^-«^ + 7^«^ + 2/3V^^) 

^yy- c/i^^^ a2(i_^2) ■ ^'-'^l 

In the no-motion limit /5 = or in the stiff fluid limit (5 = 5^ = 0), the coefficients Tr^x and 
Tyy vanish as one could have expected from the previous work |^. 
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